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ABSTRACT 

We present the discovery with WISE of a significant infrared excess associated 
with the eclipsing post-common envelope binary SDSSJ 030308.35+005443.7, the 
first excess discovered around a non-interacting white dwarf+main sequence M 
dwarf binary. The spectral energy distribution of the white dwarf+M dwarf 
companion shows significant excess longwards of 3/im. A T e ff of 8940 K for the 
white dwarf is consistent with a cooling age >2 Gyr, implying that the excess 
may be due to a recently formed circumbinary dust disk of material that extends 
from the tidal truncation radius of the binary at 1.96 R & out to <0.8 AU, with 
a total mass of ~10 20 g. We also construct WISE and follow-up ground-based 
near-infrared light curves of the system, and find variability in the K band that 
appears to be in phase with ellipsoidal variations observed in the visible. The 
presence of dust might be due to a) material being generated by the destruction 
of small rocky bodies that are being perturbed by an unseen planetary system 
or b) dust condensing from the companion's wind. The high inclination of this 
system, and the presence of dust, make it an attractive target for M dwarf transit 
surveys and long term photometric monitoring. 

Subject headings: circumstellar matter-binary stars-white dwarfs 
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Introduction 



The WISE InfraRed Excesses around Degenerates (WIRED) survey is designed to de- 
tect infrared (IR) excesses around white dwarfs (WDs) usi ng photometry from the WISE 
mission (for a complete description of the WISE mission see lWright et al.l 120101 ) . Dust, low 
mass companions, and cyclotron radiation from accreting magnetic WDs all emit at mid-IR 
wavelengths, providing a rich variety of sources to be discovered. Because of the all-sky cov- 
erage of WISE, WIRED provides a more systematic and less biased search for IR excesses 
around WDs than those performed with targeted Spitzer observations. 



In iDebes et al.l ( 1201 lal . Paper I), we charac terized the IR exce ss discovered around 
GALEX 1931+0117 using WISE photometry. In bebes et all fhoilbl . Paper II), we cross- 
correlated visible, near-IR , and WISE ph otometry of WDs discovered in the SDSS DR7 
Preliminary WD catalog (jKleinmanl |2010J) to discover excesses around WDs due to dust 
disks as well as low mass stellar and substellar companions. In the process, we have dis- 
covered a WD+dM system,SDSS J030308.35+005444.1 (hereafter SDSS J0303+0054), that 
shows an infrared excess above the photosphere of the M star companion. SDSS J0303+0054 
has a Wl — W3 color of 1.8±0.1, much redder than the average Wl — 1^3=0.5+0.2 of 
M4-M6 dwarfs detected by WISE with W3 magnitudes having signal-to- noise (S/N)> 5 
(iKirkpatrick et al.lboilh . 



SDSS J0303+0054 is an eclipsing post- co mmon envelope binary first discovered as a 
WD+dM sys t em in the SDSS DR4 WD catalog jEisenstein et al.ll2006[ Isilvestri et aDl2006h . 
Pyrzas et al.l ( 120091 . hereafter P09) detected significant radial velocity variations in the SDSS 
J0303+0054 system, as well as eclipses of the WD from optical photometric monitoring. 
They determined that SDSS J0303+0054 is a DC+dM4 binary with a 3.2 hr orbit. From the 
combination of radial velocities and light curve modeling they constrained the mass of the 
WD, assumed to have a helium dominated atmosphere, to 0.88-0.95 M and the companion's 
mass to between 0.22-0.28 M . The T e ff of the WD is not well constrained due to a lack 
of absorption features that might otherwise be detected if not for the overwhelming flux of 
the M dwarf companion; however, P09 estimated that it must be <8000 K based on spectral 
decomposition of the two stellar components. 



Observations at u, g, and i were also conducted by iParsons et al.l (120101 ). showing al- 
most complete eclipses of the WD at short wavelengths, a brigh tening of the system d uring 
eclipse of the WD at u and g, and ellipsoidal variations in i. iTappert et al.l ( I2011al ) also 
found evidence for multiple Ha components, which could be caused by a combination of 
activity /irradiation of the M dwarf and the accretion of material on the WD surface. 



Based on current prescriptions of orbital migration after common envelope (CE) evo- 
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lution (e.g. ISchreiber & Gansickd 12003 ) . SDSS J0303+0054 should have left the CE with 
a period of ~0-2 day , and slowly evolved through magnetic braking to its current orbit 



( IZorotovic et al.ll2010l ). where the M dwarf does not overfill its Roche lobe and hence no 
mass transfer through the LI point should occur in this system. Once the companion fills its 
Roche lobe and mass transfer through the LI point begins, the binary becomes a cataclysmic 
variable (CV) and circumbinary dust is a distinct possibility. Some C Vs have been found to 
show evidence for cold dust in a circumbinary disk ( Hoard et al.ll2009l ). This dust could orig- 
inate from condensation of gaseous or dusty material that is lost from the inner binary; for 
example, during classical nova outbursts, in a wind from the accretion disk and/or secondary 
star, or as "spillage" from the mass transfe r process between the two stars. Numerous single 



WDs show IR excesses due to dust (e.g. iFarihi et al.l |2010| ; Ku fc Jural l2011t iDebes et al. 



201 lbl . and references therein), likely produced by the tidal de struction of asteroids per- 



turbed into the WD Roche limit by remnant pla netary systems (IDebes fc Sigurdssonl 12002 



Jurall2003l ; iBonsor et al.ll2011t IDebes et al.ll201 id ) . Given the large mass and relatively large 
cooling age of the WD in SDSS J0303+0054 (~3 Gyr, P09, the presence of an infrared excess 
is an intriguing development in the study of this binary system. 

We present follow-up observations of the SDSS J0303+0054 system in SJ2J In Sj3]we con- 
struct a spectral energy distribution (SED) of the WD+dM in order to accurately determine 
the SED of the IR excess. We construct some simple models of the excess in £JU In §5] we 
present our follow-up near-IR observations of SDSS J0303+0054 as well as timeseries WISE 
photometry to determine the possible photometric variability in this binary, and discuss the 
implications of our work in §H1 



Observations 



In addition to GALEX, SDSS DR7, 2 MASS, UKIDSS, and WISE photometry, follow-up 
Near-IR observations of SDSS J0303-0054 were obtained in clear conditions on 20 and 21 
October 2011 with the 4.2 m William Her schel Telescope using the Long-Slit Intermediate 
Resolution Infrared Spectrograph (LIRIS; iManchado et al.lll998l ). The follow-up data were 
obtained to search for photometric variability on orbital timescales and to obtain additional 
contemporaneous photometry of the system. Additionally, we obtained the All Sky Single 
Exposure photometry^ and WISE 3-Band Cryo Single Exposure photometry^ in Wl and 
W2 of SDSS J0303+0054 in the WISE All-Sky and 3-Band Cryo Catalogs in order to track 



1 http:/ /wise2. ipac.caltech.edu/docs/release/prelim/expsup/secl_2. html#singlexp_images 
2 http:/ /wise2. ipac.caltech.edu/docs/release/allsky/expsup/sec7_l. html 
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possible short-term variability. 

On the first night of the LIRIS run, a single 9-point dither pattern of images was obtained 
in the J, H, and K s bands, with total exposure times of 18 s in each filter. This was followed 
immediately by 10 identical and consecutive sequences at H, then again at K s . On the second 
night, 50 such consecutive sequences were executed at K s only, spanning approximately 
1.2 hours in total elapsed t ime. On both nights, four standard star fields (ARNICA, 11 



stars total jHunt et al.lll998l ) were observed in a similar manner for photometric zero-point 
calibration. The data were reduced in the standard manner, by subtracting a median sky 
from each image in the dithered stack, flat fielding (using sky flats), then averaging and 
recombining frames. All data taken in the Ks-band filter were flux-calibrated using the 
ARNICA K-band standard star photometry, and the error introduced by this should be 
significantly smaller than the 5% absolute calibration uncertainty. 

LIRIS suffers from what is known as a detector reset anomaly, which appears in certain 
frames as a discontinuous jump (in dark current) between the upper and the lower two quad- 
rants. To remove this unwanted signal, after flat fielding and sky subtraction, the detector 
rows were collapsed into a median column (with real sources rejected), and subsequently sub- 
tracted from the entire two dimensional image. The resulting fully reduced frames exhibit 
smooth backgrounds, free of the anomalous gradient. 



3. Construction of the system SED 

Figure d] shows LIRIS K„ WISE Wl, W2, and W3 images of SDSS J0303+0054. It is 
detected with a S/N of >100, 43, 42, and 11, respectively. Investigation of the PSF in each 
WISE band does not reveal any extension: the science target in each band appears to be 
consistent with the expected PSF FWHM. Furthermore, inspection of higher spatial reso- 
lution UKIDSS, and LIRIS images show no indication of any neighboring sources, limiting 
contamination from comparably bright sources to angular separations <0'/35. Table [JJ lists 
the GALEX, SDSS, near-IR, and WISE photometry for this system, the modified julian date 
(MJD) of each observation after a heliocentric correction and its implied orbital phase for 
the binary, as well as our best fit model described below. 



3.1. Empirical M dwarf photometry 



An integral part of our SED modeling is the comparison to empirical SEDs of M dwarfs 
with known spectral types. In order to construct our empirical SEDs, we gathered optical, 
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NIR, and WISE photometry of known, bright M d warfs. Average op tical Sloan g — r colors 
as a function of spectral type were collected from I West et al.l (120111 ), while averag e r — i, 
i — z, and z — J colors and absolute J magnitudes were obtained from lHawley et al.l (120021 ). 
These colors were tied t o 2MASS and WISE photometry collected for M, L, and T dwarfs 
( iKirkpatrick et al.l 120111 ) . For a smooth distribution of colors to allow for interpolation 
between spectral types and to estimate the potential uncertainty in the SED templates, we fit 
low order polynomials to each color and used those fits for each color. For the u and GALEX 
filters, we assumed a Wien extrapolation from the g filter where F(X) oc (A/0.477 /zm) 5 . 



3.2. Fitting models to the SED 

In order to determine the significance of the possible nature of the IR excess, one must 
accurately subtract off the emission from the two stellar component s. We first converte d 



the photometry of SDSS J0303+0054 into units of flux density as in iDebes et al.l ( 1201 lbl ). 
and then simultaneously fit empirical SEDs of M dwarfs and model SEDs of WDs to the 
observed photometry. WD SEDs were constructed from He-rich cooli ng models kindly pro- 



vided by P. Bergeron tha t include the GALEX and WISE bands (IBergeron et al.l 11995 



Holberg fc Bergeron! 120061 ). We compared our models with GALEX FUV through J pho- 
tometry of SDSS J0303+0054, determined a median scaling between the observed and model 
photometry, and minimized the xl value to determine the best fit in logg, T e ff, d and spec- 
tral type. The resulting best-fit parameters and their 98% confidence intervals are shown 
in Table [2j Our best fit model requires a WD with logg=8.5, T e //=8940 K and an M4.5 
companion at a distance of 134 pc-the expected photometry is listed in Table [TJ The inferred 
cooling age of the WD from our models is 2 Gyr. The resulting gravity and distance to the 
WD are dependent on our assumption that the companion has the same luminosity as a 
typical field M4.5 dwarf and is thus relatively uncertain given the inherent uncertainties in 
our sample of field dwarf absolute magnitudes and colors. However, the observations of P09 
quite accurately determined both WD and M dwarf radii and masses, which we use in §4] 
and list in Tables [2] and [3J Figure [2] shows a comparison of the observed photometry and 
our best fitting model. The photometry shows excess in the H band, with all three sets of 
near-IR data showing varying levels of excess above the photosphere. However, the three 
Near-IR epochs are marginally inconsistent with each other: relative to 2MASS, the LIRIS 
H photometry is inconsistent at the 4-cr level, and the UKIDSS K photometry is inconsistent 
at the 3-<t level. These are small enough to conceivably be caused by underestimating the 
photometric errors, or could signal photometric variability in the system. We address that 
further in §5j Subtracting off the model stellar SED, we find that the Wl, W2, and W3 
bands have excesses of 22, 24, and 10 a above the expected combined stellar photospheres 
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at wavelengths > 3/im. 



Models for the IR excess 



We now investigate the origin of the IR excess around this system. IR excesses around 
presumed single WD systems arise from a) cool, low luminosity companions, b) dust disks, 
or c) cyclotron emission. Integrating over the excess and the fitted secondary SED (assuming 
a Ray leigh- Jeans approximation to the flux beyond 12/xm) results in Lm/ L*=0.3. The lack 
of any evidence for strong magnetic fields from the WD sugge sts that it is unlikel y to be due 
to cyclotron emission (e.g. as in the magnetic CV EF Eri, see iHoard et al.l ( 120071 )). The fact 
that the excess is brighter than the M companion at long wavelengths rules out a tertiary 
low-luminosity companion as the cause of the excess. We instead examine if a circumbinary, 
optically thin, dust disk can qualitatively match the observed excess. The large emitting 
surface area implied by the relatively large Lm/L* requires a circumbinary, rather than a 
circum-primary or circum-secondary disk-the implied radius of such a disk is larger than 
the orbital separation of the binary and thus would not be stable around either component 
(lArtymowicz &: Lubowl 1 1994 ). In order to model the SED of the two stellar components 
and a circumbinary dust disk, we use the procedure de veloped to model I R emission from 
circumbinary dust observed around the CV V592 Cas flHoard et al.l 120091 ). The model is 
designed to predict emission from several components i n a CV, including an optically thick 



accre tion disk and an optically thin circumbinary disk (IHoard et al.ll2007l ; iBrinkworth et al. 



20071 ). In the case of SDSS J0303+0054 there is no accretion disk and so the only constituents 
are the WD, M companion, and the dust. 

The adopted fundamental parameters of the WD and M companion that we use in our 
model are given in Table O These values are consistent, to within our fitting uncertainties, 
to the mass and gravity we derived independently for the WD from the SED alone and are 
consistent with the constraints placed by P09 from spectral decomposition and light curve 
analysis. This combination of parameters gives a statistically indistinguishable xt fr° m our 
best fitting model. 

Using these assumptions, our model for the circumbinary disk is shown in Figure [3] 
compared only to our roughly contemporaneous LIRIS and WISE data. Because of the small 
number of photometric points, we do not constrain a formal best fit-given the significant 
degeneracies present in SED modeling, we seek only to obtain a qualitative match to the 
data to determine whether the circumbinary disk hypothesis is plausible. In general, our 
simple model reproduces the observed excess, starting at H and extending out to 12 /im. 
The model underpredicts the Wl point, which could be due to structure in the disk that 
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deviates from the basic structure we assumed. The inner edge of a circumbinary disk is 
defined by the innermost stable orbit before tidal forces from the binary dom inate a dust 
grain's orbit and ejects or accretes the dust grain (lArtymowicz fc Lubowlll994f ). Therefore, 
we assume a ti dal truncation radiu s of 218 Rwd (=1-96 R Q , 1.7 times the binary separation 
as described in iHoard et al.l (120091 ) ). with an inner radius temperature of 1400 K , and an 
outer temperature of 50 K (at 18500 Rwd=0.77 AU)and a disk co-planar with the binary 
system. The inferred mass of the dust is 3xl0 20 g assuming 1/im silicate grains. If micron- 
sized silicate grains were present, however, one might also expect the W3 band photometry 
to be brighter due to a strong silicate feature at ~10 fim. This could be accounted for by 
assuming a smaller outer radius to the disk, but at this time a specific composition for the 
dust is not constrained by the data, but is used to infer a rough order of magnitude to the 
mass of the dust. Furth ermore, 10 /zm silicate features are not observed in dust around CVs 
(IBrinkworth et al.l 120071 ). implying that the dust is larger than 1/zm. Assuming the dust all 
originated from a single body composed solely of silicates, the radius would be 21 km. 



5. Near-IR LIRIS and WISE Lightcurves 

Aperture photometry of all science targets and standard stars in our LIRIS observations 
was performed using an r = 3". 75 aperture radii and sky annuli of 5" — 7". 5 in size, including 
extinction corrections. The derived JHK zeropoints for both nights agree to within 0.01 mag, 
but the scatter in the 11 standard stars each night was 0.05 mag. All photometry was 
calibrated using the J, H, and i^-band system of the standard star data. Absolute fluxes for 
SDSS J0303+0054 were derived using the derived zeropoints for each night, and also by using 
2MASS photometry for 5 relatively bright comparison stars in the field. The photometry 
thus derived for SDSS J0303+0054 using both methods agrees very well, within 0.02 mag. 

The comparison stars were also used to search for variability in time-series relative 
photometry in the K band. The flux ratio of SDSS J0303+0054 relative to the comparison 
stars was calculated for all photometry resulting from each imaging sequence, with a variance 
weighted-average of each set of resulting magnitudes for SDSS J0303+0054 to create a final 
light curve of the system. The heliocentric julian date was calculated from the midpoint 
of each photometric observation and converted to a phase using the ephemeris published in 
P09. Based on the uncertainty in the period, the phases we calculate should be accurate to 
within <10~ 2 or about 83 seconds. 

WISE observed SDSS J0303+0054 over the course of 1 day on 29 January 2010 and 
another day on 06 August 2010 for a total of 23 visits. Each visit had a detection of the 
system, with S/N of 30 and 10 in the Wl and W2 bands, respectively. We phased the data 



-8 - 



in the same manner as our i^-band data. 

Figure H shows the resulting K-, W1-, and W2-band light curves , for SDSS J0303+0054. 
In both plots, Am is negative when the system is brighter and positive when it is dimmer. 
We find that the lightcurve in the Wl and W2 photometry shows tentative evidence of 
variability, especially in a minor linear trend at phase~0.2, but the scatter in the data is 
large compared to this possible trend. The K band lightcurve is clearly variable: the system 
varies at the level of 8% peak-to-trough. Based on our models, we expect the M dwarf to 
contribute close to 82% of the total flux at H and K and 50% (41%) of the flux at Wl (W2). 



Parsons et al.l ( 120101 ) found similar variations in i-band lightcurves with an amplitude 
of 10% peak-to-trough and interpreted them to be due to ellipsoidal variations due to tidal 
deformation of the M dwarf by the WD. In z-band, the M dwarf should contribute 80% of the 
flux in the system, comparable to H and K . If both H and K vary by «10% during an orbital 
period, this could account for some of the variation observed in the near-IR photometry. The 
phase of our if -band observations overlaps with the minimum we observe at K, potentially 
explaining the dimmer value compared to 2MASS and UKIDSS. We have also overplotted 
the calculated phases for the 2MASS and UKIDSS K photometric data. The 2MASS K s 
photometry falls roughly at a phase in the orbit where the M dwarf is expected to be close to 
an inflection point of the lightcurve and close to its median brightness. If we extrapolate the 
behavior of the LIRIS lightcurve to the orbital phase of the UKIDSS photometry, we would 
expect the UKIDSS K magnitude to be ~13.31, almost 0.2 mag brighter than what was 
observed in 2005. If the UKIDSS photometry is truly dimmer, this would imply variability 
in the amount of dust in the system (or of the M-dwarf) on a timescale much longer than 
the orbital period. Observations of a full orbital phase in the if-band, as well as long term 
monitoring of this system, are necessary to confirm these unresolved issues. 

The variability at K is in phase with the ellipsoidal variations observed at i and con- 
sistent to what one would expect based on the tidal distortion of the M star from the white 
dwarf. Based on Roche geometry calculations neglecting limb darkening and using the mass 
and radius values we adopted in our circumbinary dust model, we would expect ellipsoidal 
variations peak-to-trough of 12% (which scales to 10% if the secondary contributes 84% of 
the flux at K). The shape of the lightcurve in Figure H] deviates slightly fro m a pure sine 



function which is also seen at shorter wavelengths (P09 iParsons et al.l 120101 ) . We further 



discuss our results in the context of all the lightcurves for this system in 
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6. Discussion 

The WISE observations of SDSS J0303+0054 represent the first likely detection of dust 
in orbit around a non-interacting post-common envelope binary. If the dust disk is aligned 
with the orbit of the binary, SDSS J0303+0054 could be a fruitful target for high cadence 
near-IR transit searches. Our photometry was sensitive to transits of depths > 50 mmag 
and durations >3 minutes. With the same photometric accuracy and sampling at < 30 sec, 
planets with radii of ~5 R e could be detected if they orbited at the inner edge of the disk 
and transited with durations of ~1 minute. 

The presence of ellipsoidal variations deviating from pure sinusoidal behavior, now seen 
at K with our new observations, also opens up the possibility that the structure of the cir- 
cumbinary disk can be directly probed by detailed analysis of long wavelength time series 
photometry. The circumbinary disk is not symmetrically heated throughout an orbital pe- 
riod, and if there exist significant azimuthal density or compositional variations, these could 
be detected as asymmetries in the lightcurve. Since the disk contributes 15% of the light at 
K, such variations could be detected at the few percent level. 

One possibility is that SDSS J0303+0054 is the binary equivalent of isolated WD dusty 
disks, which most likely form from the tidal disruption of rocky bodies a few tens of kilometers 
in radius. Given the relatively long cooling age of this system compared to other dusty WDs 
and the inferred T e ff of the WD, this would be comparable to the coldest recorded dusty 



single WD observed, G166-58 (IFarihi et al.ll2010[ ). The amount and location of dust present 



at SDSS J0303+0054 would not be detectable around a single cool WD. The presence of 
the M companion's additional luminosity, which far exceeds that of the WD, enhances the 
detectability of the circumbinary dust. 

If the dust is caused by a tidally disrupted asteroid, then SDSS J0303+0054 is reminis- 
cent of Kepler-16b, which contains a 0.6 M m primary, a 0.2 M & companion, and a Saturn- 



mass companion in a circumbinary orbit (IDoyle et al.ll201l[ ). Such a binary, with a 41 day 



orbit, would eventually evolve into a common envelope binary. The remains of the circumbi- 
nary planetary system could feed a similar dust disk to that seen around SDSS J0303+0054. 
However, several potential problems exist for such a system. Firstly, the Poynting- Robert son 
drag timescale mi ght be prohibitively short for such a disk. The P-R drag timescale can be 



approximated by ([Hansen et al.l 120061 ): 



Ps • > ■ ■ • - 



TpR - 10 ° ( >m7 Ug~^J llQiW (lF^ 1 ^ 



where is the luminosity of the star, s and p s are the average grain size and density 
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respectively, and r is the distance from the central star. As noted previously, the M dwarf 
dominates the luminosity of the system. Assuming a bolometric luminosity of 5xlO _3 L Q 



(IRebassa-Mansergas et al.l 120071 ). and assuming a minimum distance of the disk from the 
M dwarf of 0.76 R & the PR timescale can have a range of values since the dust orbits the 
binary's center of mass. At the closest distance to the M dwarf, the P-R timescale is 6 yr, 
while at the furthest distance away from the dust disk inner edge the P-R timescale can be as 
long as 100 yr. The relatively short P-R timescale implies that dust must be replenished at a 
rate of 5xl0 19 g/yr or 2xl0 12 g s -1 . The minimum replenishment rate, corresponding to the 
longer P-R timescale, would be sixteen times smaller. This rate of dust production is high 
compared to the inferred accretion rates of dust onto single WDs, but could be reconciled if 
the radius of a typical dust particle is larger than 1/im. 

The mechanism for creating the dust in the asteroidal scenario might be significantly 
different than for single WDs. It is unclear how an asteroid could tidally disrupt at the 
circumbinary disk inner radius, which is much larger than the Roche disruption radii of 
either the WD or the M dwarf. If the dust is from tidally disrupted bodies, multiple asteroids 
could have similar periastra, and the dust forms from multiple collisions of these objects. 
The accretion rate would be lower if the disk were optically thick, since P-R drag could not 
act on all of the dust present. Assuming an inclination of 82°, an inclination consistent with 
the results of P09, the lack of reddening toward the system constrains the vertical height of 
the disk to be < 0.1 R . 

Alternatively, low levels of mass transfer that might be present in this system could 
account for the presence of a dust disk as well. The combined photometric and spectroscopic 
behavior of the SDSS J0303+0054 system suggests low levels of wind material being accreted 
by the WD. The u and g light c urves of the system show a brightening of the system 



coincident with the WD's eclipse. iParsons et al.l (120101 ) interpreted this excess flux as due 



to the side of the M-dwarf facing away from the WD being brighter. However, this is 
contradicted by the ~96% depth of eclipse seen in the u band, despite a 10-20% increase 
in flux leading up to the eclipse. If the M companion were responsible for the excess, the 
eclipse would be shallower than that predicted by the ratio of the M dwarf optical flux to 
the total flux from the system. Based on our best fit model, the eclipse depth should be 
~90% in u. This difference suggests our model over-predicts the flux from the M dwarf, or 
under-predicts the flux from the WD. We instead interpret the brightening, which covers 
about half the orbit, as due to the side of the WD facing the secondary being hotter due to 
broad wind accretion on t hat hemisphere. The presence of multi-component Ha emission 



lines (ITappert et al.l 1201 lai ) are also indicative of accretion from the secondary wind. 



Similar to CVs with infrared excesses, this wind accretion may provide enough "spillage r 
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to create a low mass dust disk. Using the above P-R drag mass loss rate as an estimate of dust 
destruction, we can infer the rough mass injection required to obtain a steady state dusty 
disk. If we assume that all of the metals in the secondary's wind condense into dust, and 
that half of the wind must contribute to the "spillage" , this suggests that the M dwarf must 
be losing mass at a rate of 5xl0 21 g yr _1 with a minimum rate an order of magnitude lower 
for the longer P-R drag and assuming roughly solar metallicity. This could be uncertain by a 
factor of 10 if the disk is radially optically thick-this can slow the destruction of dust grains 
if some of the dust is shadowe d. If we com pare this to the mass loss rate of other M dwarf 
winds that accrete onto WDs JPebesll2006L we find that the M dwarf in SDSS J0303+0054 
must be losing mass at a maximu m rate ~100-1000 times that of typical M dwarfs and m ore 



than the solar wind mass loss rate (jPebesI l2006t iTappert et al.ll2011bl ; iPyrzas et al.ll2012l ). If 
the M dwarf is tidally distorted, as required by the presence of ellipsoidal variations, then 
this could enhance the mass loss in its wind. The escape velocity at the secondary's surface 
is smaller than the escape speed of the binary system, which could naturally trap a majority 
of the wind particles in stable orbits around the binary. 

Neither scenario completely explains the origins of this unusual system. If there is sig- 
nificant dust production from low levels of mass transfer or from remnant planetary systems, 
it is possible that a significant number of other PCEBs contain dust as well. With over 1000 
candidate WD+M systems from the WIRED survey, a more careful analysis of that popu- 
lation may yield other dusty PCEBs, and hopefully lead to a resolution to the origin of the 
dust, and how it fits into the evolution of compact binaries. 



This work is based on data obtained from: (a) the Wide-field Infrared Survey Ex- 
plorer, which is a joint project of the University of California, Los Angeles, and the Jet 
Propulsion Laboratory (JPL), California Institute of Technology (Caltech), funded by the 
National Aeronautics and Space Administration (NASA); (b) the Two Micron All Sky Survey 
(2MASS), a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center (IPAC) /Caltech, funded by NASA and the National Science Foundation; 
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Table I. SDSS J0303+0054 Photometry 
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a Where multiple observations exist, a mean value is calculated 
b Phase calculated from ephemeris of P09 
C AB magnitude, rather than Vega magnitude 
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Table 2. SDSS J0303+0054 Fitted Fundamental Parameters 



Fundamental Parameter Value P09 Constraints 



WD log g 
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WD Mass (M ) 
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Inferred Distance (pc) 
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Table 3. SDSS J0303+0054 Adopted Fundamental Parameters 



Fundamental Parameter White Dwarf M dwarf 



Mass (M ) 


0.89 


0.25 


Radius (R©) 
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Inferred T e ff 
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Inferred Distance (pc) 
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Luminosity (L ) 
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Fig. 1.— Infrared images of SDSS J0303+0054 in LIRIS i^-band, and the WISE Wl, W2, 
and W3 bands centered on its SDSS DR 7 coordinates. North is up and East is to the left. 
Greyscale goes from 2a below local background level to 25a above the local background. 
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Fig. 2.— SED of the SDSS J0303+0054 system from GALEX FUV to W3, with an upper 
limit in W4. Black squares are the observed photometry, including the near-IR photometry 
at J, H, andK at three separate epochs. Errors are typically smaller than the symbol size. 
Overplotted with red asterisks is a best-fitting model DC+dM6 SED, with a T e ^=8940 K, 
log g=8.5 WD. 
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Fig. 3. — Comparison of SDSS J0303+0054's photometry to our model circumbinary dust 
disk. The black line is the combination of the WD (blue line), M4.5 companion (orange line), 
and circumbinary dust disk (red line) consisting of 1/xm grains arranged in a geometrically 
thin, optically thin disk with an inner radius of 1.96 R & and inclination co-planar to the 
best fit inclination from P09 (80°-84°). The outer radius of the disk is not well constrained 
by our models, but may extend as far as 0.8 AU. The M dwarf companion has a luminosity 
an order of magnitude greater than the WD. 
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Fig. 4. — (Top) K light curve (A m=m-13.34) of SDSS J0303+0054 phased to the ephemeris 
of P09. Approximate uncertainties in the phasing should systematically be ±0.01. The dia- 
mond overplotted at phase=0.8 is the epoch 2005 UKIDSS K photometry, and the diamond 
at phase=0.17, the epoch 2000 2MASS K photometry (similar phasing uncertainties, for 
more discussion see §5]). Variability in the i^-band brightness is evident, and we overplot 
a sinusoid to the data assuming a period half that of the orbital period, an amplitude of 
37 mmag, and zeropoint at 0, consistent with the ellipsoidal variations observed at shorter 
wavelengths. This is meant to guide the eye, rather than represent a true fit to the data. 
(Bottom) WISE Wl (black asterisks) and W2 (red diamonds) photometry from 2010 and 
2011 (Am = m — m wi w2) phased to the same ephemeris. There may be tentative vari- 
ability in the data. 



